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abstract 

Interaction  between  the  oceanic  vertical  shear 
field  and  the  natural  response  of  a  recently 
developed  expendable  dissipation  profiler  (XDP) 
causes  errors  to  be  introduced  in  small  scale 
velocity  measurements.  A  numerical  model  was 
developed  to  study  this  problem  and  two  sources  of 
errors  were  determined,  the  primary  source  being 
the  shear  induced  torques  at  the  instrument  length 
scale.  ^ 


1.  INTRODUCTION 

The  direct  measurement  of  the  dissipation  rate 
of  small  scale  turbulence  has  long  been  a  goal  of 
oceanographers,  but  has  become  feasible  only  in 
recent  years  with  the  development  of  probes  capable 
of  measuring  shears  on  subcentimeter  scales.  The 
primary  sensor  for  this  type  of  measurement  is  the 
air-foil  shear  probe,  initially  developed  by  Ribner 
and  Siddon  <1969,  cited  in  Osborne  and  Crawford, 
1978)  and  subsequently  improved  for  water  operation 
by  Osborne  and  Crawford  (1978). 


The  airfoil  probe  operates  on  the  pr 
that  an  ogive,  passing  at  constant  velocity 
a  flow  field  with  shear  at  right  angles 
motion,  will  experience  a  lift  force  which 
function  of  the  shear  field  and  its  veloc 
the  ogive  is  elastic,  it  will  deform  in  a 
nearly  proportional  to  that  lift  force.  The 
probes  are  constructed  by  mounting  a 
of  piezoelectric  material  in  a  silicone 
ogive  that  is  mounted  on  the  end  of  a  long 
at  the  front  of  the  vehicle. 
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The  silicone  rubber  is  compliant  and  imparts 
any  deformation  to  the  beam,  the  electrical  output 
of  which  is  amplified  and  filtered  by  the  circuitry 
in  the  vehicle.  The  data  may  be  stored  internally 
in  the  vehicle  as  in  the  Camel  (Osborn,  1977);  or 
they  may  be  sent  to  the  surface  via  a  direct  data 
link  as  in  the  Advanced  H ic ros t rue t ure  Profiler 
(AMP)  (Gregg,  1982)  or  the  expendable  dissipation 
profiler  (XDP). 

Tne  XDP  was  designed  and  built  by  Lueck  and 
Osborn  (1983)  as  a  relatively  inexpensive 
expendable  probe  to  measure  dissipation  from  a 
moving  vessel.  The  basic  design  of  the  vehicle  is 
similar  to  the  well-known  XBT,  with  certain 
differences  that  were  int.-oduced  to  reduce  the 
mechanical  vibration  of  tiie  instrument. 


One  phase  of  the  Ocean  Measurement  Sensors 
Project  at  NORDA  has  been  the  evaluation  of  XDPs. 
In  November,  1984,  12  XDPs  were  deployed  during  the 
TROPIC  HEAT  cruise  of  the  University  of  Washington 
Applied  Physics  Laboratory.  Each  of  the  XDPs  was 
dropped  immediately  after  UU/APL  had  completed  an 
AMP  (Advanced  Microstructure  Profiler)  drop. 

The  AMP  (Gregg,  1982)  has  two  airfoil 
turbulence  probes  in  its  sensor  suite.  The  AMP 
vehicle  is  longer  and  more  massive  than  the  XDP 
vehicle  and  is  expected  to  provide  a  much  more 
stable  platform  from  which  to  measure  turbulent 
velocity  fluctuations.  This  is  confirmed  by  the 
excellent  agreement  with  the  Nasmyth  spectral  form 
(Oakey,  1983)  for  the  observed  level  of  dissipation 
(obtained  from  integrating  the  shear  spectrum). 

When  the  records  from  both  instruments  were 
compared,  the  XDP  showed  much  larger  apparent 
dissipation  rates  than  those  obtained  from  the  AMP. 
The  XDP  spectra  often  were  found  to  have  unexpected 
peaks  in  1  to  10  cpm  wavenumbei  range  and  an  excess 
of  variance  at  high  wavenumbers  over  the  Nasmyth 
spectra  for  the  same  integrated  shear  variance.  A 
typical  spectrum  of  the  .XDP  shear  data  is 
illustrated  in  figure  1. 

If  the  XDP  is  to  be  used  as  a  reliable 
in.strument  for  turbulence  measurements.  the  cause 
of  the  differences  between  tiie  AMP  and  XDP 
measurements  must  be  located  and  corrected.  Because 
some  of  the  peaks  occurred  in  the  1  to  10  cpm  band, 
we  suspected  that  the  background  shear  field  might 
be  amplified  by  the  natural  modes  of  oscillation  of 
the  XDP  vehicle  at  scales  near  the  vehicle  length. 

It  is  difficult  and  expensive  to  design  and 
conduct  an  experiment  to  measure  the  amplification 
of  the  actual  shear  field  in  the  ocean  by  an  XDP.  A 
numerical  simulation  offers  a  simpler  means  of 
testing  this  hypothesis  and  piovides  additional 
information  that  could  not  be  obtained  as  regards 
the  mechanics  of  the  XDP  vehicle  without  a  major 
engineering  effort. 


2.  METHOD 

Tfie  basic  method  followed  ii,  peiforning  this 
simulation  is  to  model  the  motion  of  an  XDP  vehicle 
by  applying  Newton's  thiid  law  in  two  degiees  of 
tianslaiion  atid  one  of  loiatM'i;  to  the  XDi'  tiody. 
The  najoi  problem  in  simulating  the  motion  of  the 
veltirle  is  to  obtain  a  good  i  opt  esen  t  a  t  i  on  of  the 
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Figure  1.  Spectru”  of  XDP  shear  data  obtained 
during  the  Tropic  Heat  cruise  of  November,  1984. 
The  data  used  for  tu.s  spectrum  vere  from  record 
XOrTHOlC.  data  poitus  1  ibO'N  using  Hanning 
-indoving  and  50.'^  overlap.  The  '^issip^i'ion  laie  vas 
computed  by  integrating  to  the  cutoff  vavonumbei . 
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ob'’ained  by  itetation.  Th^.-  superposed  curve  is  tf»e 
Nasmyth  spectru.m  for  rh^  estim.arod  dissipation 
1  a  t  e . 

forcing  field:  the  detailed  velocity  stiuctuie  n: 
•he  ocean  as  a  fur.'’ lion  of  depth. 

In  this  model,  the  vehicle  is  constrained  to 
planar  motion.  This  assumption  vas  made  both  fot 
computational  simplicity  and  ease  of  interpretation 
of  the  results.  The  othe:  main  assumption  vas  the 
value  of  the  virtual  ma  and  moments  of  inertia. 
The  virtual  mass  vas  taken  a'-  the  sum  of  ilie 
vehicle  mass  and  that  oi  the  vater  displaced  by  the 
vehicle.  The  virtu.jl  moment  of  inertia  vas  taken  as 
the  sum  of  the  vehicle  moment  of  inertia  about  the 
center  of  gia'-ity  and  Mie  moment  of  inertia  of  the 
displaced  water.  The  coordinate  sy.stem  for  the 
instrument  model  is  shown  in  figure  2. 
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fiRuip  7.  Coor  iinate  system  for  XOF  vehicle.  The 
.ifiKle  a  is  t.vKen  positive  in  the  counterclockwise 
vense.  Tlie  orip.in  foi  {,  is  the  center  of  gravity 
tT.ass:)  of  the  XDP  vehicle. 

:  c.  t  a  1 1  on  ) .  .and  .1)  the  inside -out  ‘-inding  of  the 

V  lie  c;  o  j  )  . 

Tfie  physical  dimensions  and  properties  of  t'n 
veliiclc  ate  summarized  in  table  1,  along  with  othei 
'  niictants  used  in  the  numerical  simulation. 

The  mass  of  the  vehicle  vas  determined  by 
veighing  the  XDP  on  a  laboratory  balance.  The 
moment  of  inertia  uas  determined  by  mounting  the 
XDP  as  a  compound  pendulum  and  measuring  the 
pet iod  (Timoshenko  and  Young,  1940,  p  401).  The 
tnf>.meni  of  inertia  about  the  centei  of  gravity  uas 
(onp.ited  by  the  usual  forruKa  (loc.  cit..  p  508). 

r..  UIFFERE.-.TI.AL  PQUAriOriC 

The  model  d  i  f  f  ereti  r  i  ,i  1  etpi.at  ions  are  derived 
from  Nevton's  tfiiid  la-.  b.ilaiu;ing  r:asr.  times 
acccleiat  ion  in  tlie  hoiizont.al  ’and  vertical 
directions  and  the  angulai  momentum  foim  of  the 
thiid  lav  for  lot.ition  in  the  y  -  d  pl.ane.  In  all 
tlii-:  vot  k ,  the  motion  of  the  instrument  is  assumed 
to  lie  const  I  a  iiied  to  oi...  plane  (v  .  II).  The 
coord i na te  system  is  cielioel  in  the  (x.r)  plane 
vitii  r  taken  poitive  upvat.i.  '-'e  .  .n,  .i-.i,,  I'.i.v 

vii!i.  ilir.  dilfeienli.il  eipiai'.  .o;  of  mr.iion  .r. 

vet  t  i  1  a  i  .  to  i  tv  . 
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a  .  ♦  •  auiguLar  velocity  , 

V  .  gravitational  acccelera t ion  ♦  buoyancy 
acceleration  +  drag  acceleration, 

u  =  horizontal  drag  acceleration, 

=  gravitational  restoring  torque  ♦ 
shear  torque 

The  specific  formulation  of  these  accelerations  is 
given  as 

-1  L  2 

w  =  I-  rag  +  b  -  pCp  (  tia  e 

LD  I  sina  |  )  v | w | | 


ii  .  ra^^CppD  L  (U{z)  -  u)|U(z)  -  uj 


^  -  bXl*  sin  a  - 


j^^tfeos  o  sin  a  0^^^)  dt, 


0^*^  -  0.5  CpPCOS  a  lU(z*J,cos  o)-u| 


|U(z*f,cos  a)-u|D(£.), 


.  -0.5  CpcDff;)  w  |w|  sin  a  . 


Vaiiahle  Description 


Angle  of  XDP  vehicle  to 
the  vertical,  positive  in 
counterclockwise  sense 

Radius  of  XDP  nose  C 

Buoyancy  1 


Transverse  drag 
coefficient 


Diameter  of  Instrument  0.05  m 


Acceleration  of  gravity  9.8  m  s 
Virtual  Moment  of 


Inertia 


Total  length  of  XDP 
Distance  CG  to  nose 


Distance  CG  to  tall 


0.0156  kg m 


0.340  m 


0.243  m 


Distance  CG  to  center  of 
buoyancy  0.073  m 


Mass  of  XDP 


1.088  kg 


Background  shear  field,  U(z) 

The  driving  force  for  this  set  of  differential 
equations  is  provided  by  the  effect  of  the 
background  horizontal  velocity  field,  U(z)  exerting 
a  drag  force  on  the  XDP  vehicle.  The  velocity  field 
was  estimated  by  constructing  a  depth  series  of 
shear  based  on  the  Garget,  et  al.  (1980)  model 
shear  spectrum.  The  specific  spectrum  used  was 
given  by: 


,  k  c  kg.  0.1  m' 


Virtual  Mass  2.176  kg 


Density  of  Seawater  1025  kg  m 

Shear  Spectrum 
Horizontal  XDP  velocity 


Background  velocity  field 
Vertical  XDP  velocity 


Cq  (k/kg) 


.  kp<  k  <  kj  .  1  m  . 


Nasmyth  Spectrum,  k  >  kj  ; 


(N/Nq)^  x2xl0^  s^/epm 


1.4  eph 


Horizontal  XDP  position 

Vertical  XDP  position 

Angular  velocity  of  XDP 

Integration  variable 
(Distance  from  center  of 
gravity,  positive  aft) 


In  the  range  k  >  1  m“\  the  Nasmyth  spectrum  The  model  shear  field  was  constructed  by 

(Oakey,  1983)  was  used  with  a  dissipation  rate  of  c  taking  the  square  root  of  the  spectrum,  multiplying 

it  by  a  complex  factor  exp|i^(k)l,  where  iCk)  is  a 
»  1  X  10  W/kg.  random  variable,  uniformly  distributed  over  (0,2fi). 
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An  inverse  Fourier  transform  was  then  applied  to 
the  resultant  function  of  k  to  get  random  phase 
shear  profile.  The  shear  profile  uas  then 
integrated  numerically  to  obtain  the  velocity 
profile. 

Symbolically  (in  continuous  form); 


U(z)  .  Jq  F'^^S^''^k)  expl  i^(k)  1  1  dC, 


tigures  3  and  U.  Figure  5  presents  the  horizontal 
motion  of  the  probe  vith  depth;  figure  6,  the 
horizontal  velocity;  figure  7,  the  rotational 
angle;  figure  8, the  angular  velocity;  and  figure  9, 
the  vertical  fall  rate.  The  "sensed"  horizontal 
velocity  was  computed  by  combining  the  velocity  due 
to  rotation  with  the  difference  in  horizontal 
motion  of  the  vehicle  and  the  background  flow; 


U 


sensed 


U(z)  -  u  -  Lj+  . 


where 

F^'{f(k)l  .  (2n)'^  exp(ikO  dk 

is  the  inverse  Fourier  transform. 

A  63336  point  discrete  Fourier  transform  was 
used  providing  a  1  cm  spacing  of  points.  As  the 
aperture  of  the  probe  was  about  2  cm,  this  was  felt 
to  be  sufficient  to  model  the  response  of  the 
probe.  For  the  rotational  torque  calculations,  this 
spacing  gave  about  34  points  for  integrating  the 
forces  over  the  length  of  the  vehicle. 

3.  DISCUSSION 

Forcing  mechanisms 

Two  mechanisms  should  affect  the  motion  of  the 
vehicle.  Both  produce  torques  on  the  vehicle  body 
yiiich  will  cause  it  to  oscillate  as  it  falls 
through  the  water  column  and  thus  induce  spurious 
additions  to  tlie  observed  transverse  current  field. 
The  first  and  most  obvious  torque  is  caused  by  the 
shear  field  nf  the  background  flow.  This  field 
will,  in  general,  induce  a  torque  because  a  varying 
distribution  of  force  will  be  applied  to  the  body 
along  its  length. 

The  source  of  the  second  torque  is  not  as 
obvious.  If  the  vehicle  were  to  fall  into  a  region 
of  constant  flow.  there  would  be  a  uniform  force 
over  the  whole  vehicle  due  to  the  difference  in  the 
horizontal  velocity  component  of  the  vehicle  and 
the  background  flow.  In  most  free  fall  vehicles, 
the  center  of  gravity  (mass)  is  below  the  center  of 
projected  transverse  area  of  the  vehicle  and,  as 
the  moment  of  force  is  computed  about  the  center  of 
gravity,  a  uniform  force  distribution  will  result 
in  a  torque.  Thus,  at  low  frequencies  (and 
wavenumbers),  we  should  expect  to  see  the  angle  of 
the  vehicle  correlated  with  the  background 
horizontal  velocity.  This  will  not  cause  problems 
with  induced  velocities,  due  to  the  low  frequency 
of  the  motions,  but  wHl  cause  an  error  due  to  the 
angular  sensitivity  oT  the  probe  (not  modeled  ir 
(his  study). 

Mo.K-1  Re's'll  ' 

The  differential  equations  were  iiuegiaied 
using  the  IMSl.  version  of  the  Dear  solver  for  stiff 
'irdiuaiy  differential  equations  (IMSL,  1982).  The 
inituil  renditions  on  .all  variables  were  set  to 
zero,  r;o  only  the  external  foicing  would  affect  the 
motion  of  the  XDP  vehicle. 

The  velocity  and  .shear  field  used  to  drive  the 
probe  .'ncjtion  .ire  plotted  as  .i  lunctinn  of  depth  in 


For  this  study,  no  attempt  has  been  made  to  model 
the  dynamic  response  of  the  probe.  There  should  be 
little  disagreement  between  the  sensed  flow  as 
computed  above  and  the  sensed  flow  after  applying 
the  probe's  low  pass  response  filter  for  the  scales 
under  considerat ion  as  the  probe  cutoff  wavenumber 
is  about  1  cm. 
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Figure  3.  Background  horizontal  velocity  versus 
depth.  The  velocity  was  obtained  by 'integrat fng  the 
model  shear  field. 


The  "spn.sed"  velocity  is  plotted  versus  depth 
in  figure  10  .  When  this  figure  is  compared  with 
figures  3  and  8,  it  appears  that  the  component 

of  the  "sensed"  velocity  is  nearly  as  important  as 
the  U(z)  component.  The  scale  of  the  plot  makes  it 
difficult  to  See  the  details  of  the  vehicle  motion 
effects  on  the  "sensed"  velocity. 

To  determine  the  coa t r Ibut ions  of  the  vehicle 
oscillations  on  the  observed  velocities  and  the 
time/space  scales  at  which  effects  Occur,  spectra 
of  both  the  "sensed"  velocity  and  the  background 
velocity  were  computed.  These  spectra  were  computed 
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Figure  12.  Transfer  function  betveen  U(:)  and 

(2).  The  sensed  velocities  are  amplified  with 

sed 

respect  to  the  background  velocities  in  the  to 

5  cpm  wavenumber  band.  The  maximum  amplification 
occurs  at  about  1.3  cpm  vith  a  30  fold 
ampl  i  f  i  c.a  t  i  on  . 


for  times  and  depths  after  which  the  probe  had 
reached  terminal  velocity  (about  1  to  3  meters). 
These  spectra  are  plotted  in  figure  II.  As  can  ho 
seen,  the  "sensed*'  velocity  field  has  a  very  .strong 
peak  in  the  I  to  10  cpm  band.  indicating  a  vety 
strong  amplification  of  the  background  shear,  as 
was  hypothesi2ed.  The  transfer  runction  amplitude 
which  is  equal  to  the  ratio  of  the  output  spectrum 
to  the  input  spectrum  was  also  computed  and  is 
depicted  in  figure  12.  This  figure  shows  that  the 
background  field  is  amplified  by  about  a  factor  of 
5  -  10  in  the  1  to  10  cpm  band. 

The  numerical  simulation  was  repe.ited  loi 
several  different  values  of  E  ranging  tron  10 

U/kg  to  U/kg  to  determine  the  effect  nf 

varying  dissipation  rates  on  the  dcgiee  of 

amplification.  These  effects  indicate  a  ■-■cak  non 
linearity:  the  amount  of  amplification  inrieases 
weakly  with  the  dissipation  ijte. 

Unresolved  Problems 

While  thi.s  numerical  study  has  .shown  rh.it  the 
probe's  tcsponse  significantly  affect::  ihi»  observed 
shear  in  the  1  10  cpm  band,  it  has  not  rlaiified 

the  I easons  for  the  large  apparent  differences  iu 
observed  i  ss  i  pa  t  io,-  rates  obtained  f  n>m  the  Xf>C 
and  the  AMP.  rii  .  A.U.  Oteeri.  NDPOA.  snnM.i  1 


communication,  has  suggested  that  some  of  the 
observed  peaks  in  the  low  frequency  band  (1  -  10 
cpm)  may  be  due  to  slow  rotations  of  the  vehicle 
coupled  with  precessional  modes  of  oscillation.  The 
high  frequency  flattening  of  the  spectrum  may  be 
due,  in  part,  to  despooling  noise  associated  with 
the  signal  wire  and  to  wake  noise  transmitted 
through  the  instrument  to  the  sensor  element. 

6.  CONCLUSIONS 

While  several  observations  of  the  XDP  probe's 
behavior  that  are  still  unaccounted  for,  this 
modeling  study  has  shown  that  there  may  be 
considerable  amplification  of  the  background  shear 
field  in  the  wavenumber  band  between  1  and  10  cpm. 
This  is  the  same  band  where  the  peak  in  turbulent 
energy  dissipation  is  expected  to  occur  (based  on 
the  Nasmyth  spectrum).  Neglect  of  the  probe 
response  in  this  region  will  cause  serious  errors 
in  estimation  of  the  dissipation  rate.  Correction 
of  the  amplification  effect  is  made  more  difficult 
by  the  nonlinearity  of  the  effect,  the  occurrence 
of  spectral  peaks  not  predicted  by  the  model  and 
high  frequency  flattening  of  the  spectrum.  It 
appears  that  further  work  is  required  before  the 
expendable  dissipation  profiler  is  ready  to  be  used 
as  a  routine  survey  tool. 
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